Nonlinear modes in binary bosonic condensates with the pseudo-spin-orbital coupling 
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We consider a binary Bose-Einstein condensate (BEC) witli nonlinear repulsive interactions and 
linear spin-orbit (SO) and Zeeman-splitting couplings. In the presence of the trapping harmonic- 
oscillator (HO) potential, we report the existence of even, odd, and asymmetric spatial modes. 
They feature alternating domains with opposite directions of the pseudo-spin, i.e., anti- ferromagnetic 
structures, which is explained by the interplay of the linear couplings, HO confinement, and repulsive 
self-interaction. The number of the domains is determined by the strength of the SO coupling. The 
modes are constructed analytically in the weakly nonlinear system. The dynamical stability of the 
modes is investigated by means of the Bogoliubov-de Gennes equations and direct simulations. A 
notable result is that the multi-domain-wall (DW) structures are stable, alternating between odd 
and even shapes, while the simplest single-DW structure is unstable. Thus, the system features 
a transition to the complex ground states under the action of the SO coupling. The addition of 
the Zeeman splitting transforms the odd modes into asymmetric ones via spontaneous symmetry 
breaking. The results suggest possibilities for switching the binary system between states with 
opposite (pseudo) magnetization by external fields, and realization of similar stable states and 
dynamical effects in solid-state and nonlinear-optical settings emulated by the SO-coupled BECs. 

PACS numbers: 67.85.-d, 03.75.Kk, 03.75.Mn, 05.30.Jp 



I. INTRODUCTION 

An important application of Bose-Einstein condensates 
(BECs) is their use for emulating a variety of fundamen- 
tal effects originating in the realm of condensed-matter 
physics, in a form which is much easier to handle in 
atomic gases [l|. Much attention was recently attracted 
to the implementation of linear couplings between two 
components of a binary BEC, which emulate the spin- 
orbit (SO) interactions in solid-state settings, with the 
respective spinor order parameter mapped into the two- 
component wave functions of mixtures of different states 
of the same atomic species • Besides the fundamental 
interest concerning the dynamics of spinor BECs, the sys- 
tem represents a tcstbed for the study of artificial gauge 
fields, which is another topic of a rapidly growing inter- 
est Q. To a great extent, the possibility of emulating 
condensed-matter effects in ultra-cold gases was brought 
to the focus of the current research by experiments Q 
demonstrating the action of Abelian and non-Abclian 
synthetic gauge fields in BEC [H. 

In addition to simulating condensed-matter phenom- 
ena, the studies of SO-coupled BECs reveal new matter- 
wave effects, produced by the interplay of the SO cou- 
pling and the mean-field nonlinearity induced by atomic 
collisions. Such effects include sophisticated vortical Jg 
and monopole structures, multi-domain patterns U, 
and patterns produced by long-range interactions J9| , tri- 
critical points [l3l, skyrmions flli. solitons [13: Ell: in- 
teraction with optical lattices [Tj], etc. Note that the 
system of coupled Gross-Pitaevskii equations (GPEs) de- 



rived in Ref. [i2[ is tantamount to that describing the 
co-propagation of two polarizations of light in twisted 
nonlinear optical fibers [15| , hence an additional link be- 
tween completely different physical settings is the pos- 
sibility to emulate birefringent optical fibers by the SO- 
coupled condensates, and vice versa. Furthermore, a re- 
cent report on the optical implementation of a model sim- 
ulating massive Dirac fermions [l^ indicates the further 
relevance of results reported below to nonlinear guided- 
wave optics. 

The variety of expected effects in SO-coupled BEC may 
be greatly enriched when three (or more) atomic levels 
are employed to define relevant components and inter- 
actions between them In particular, tripod atomic 
schemes give rise to coupled GPEs featuring various non- 
linear terms which do not amount to the usual self-phase 
and cross-phase modulation [l3l . Such systems open pos- 
sibilities for nonlinear control of SO-coupled BECs by 
means of resonant laser illumination. 

While spinor BECs typically give rise to two phases 
with opposite polarizations of pseudo-spins, and to mix- 
tures of such phases 0, H, not any form of the SO 
coupling results in energy splitting between the phases. 
For such a situation, the role of the nonlinearity is cru- 
cially important, as it lifts the degeneracy and induces 
the phase separation, cf. Ref. [T8| . 

In this work we study nonlinear modes originating from 
degenerate linear eigenstates in the SO-coupled binary 
BEC loaded into a harmonic-oscillator (HO) trap, finding 
a full set of such nonlinear states. We produce basic spa- 
tial patterns which exist in the case when the Rashba [l^ 
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and Dresselhaus [20[ couplings have equal strengths. The 
patterns are built of alternating domains with opposite 
directions of the pseudo-spin, with each state character- 
ized by zero or nonzero total (pseudo-) magnetization. 
The number of domains is determined by the strength of 
the SO coupling. The most essential new results concern 
the stability of the competing states, which is determined 
by the total magnetization. In contrast to previously 
studied settings, we demonstrate that the multi-domain- 
wall (DW) patterns may be stable while the single DW is 
not, hence the system's ground state shifts to the complex 
patterns. We also report a possibility of the dynamical 
switching between different stable patterns. 



II. THE MODEL 

We consider the two-component effectively one- 
dimensional (ID) BEG described by spinor ^ ~ 
col{'ipi,'ip2}, whose components ipi and "02 represent 
pseudo-spin components | t) and | -I). The dynamics 
of the system is governed by the Hamiltonian, H = 
Hq + Hirit, where, in scaled units, Hq ~ J^!^ '^^'H'^dx, 
n = (1/2) (-92 + + flas) + inaida;, with Pauh ma- 



trices (73,1, and X /2 is the axial trapping potential [21 



Zeeman splitting f2 is induced via a constant magnetic 
field acting along the z-axis, while the SO coupling, ac- 
counted by coefficient k, results from a combined effect of 
the Rashba and Dresselhaus couplings, and is determined 
by intensities and wavelengths of laser beams which cou- 
ple the relevant atomic levels. The interaction terms 
for the underlying energy-level scheme are represented 
by iJi„t = /+r[(5i/2) + |V;2|^) +5|^inV2p] dx 

p7| . The derivation of the ID model from the full 3D 
one follows the scenario of introducing a tight HO po- 
tential of the transverse confinement and factorizing the 
wave functions into the ground state of that potential and 
a free longitudinal wave function, which has been elab- 
orated in detail [g^]. The structure of the SO-coupled 
system does not imply any problem in following this sce- 
nario, provided that it leads to the usual cubic nonlinear- 
ity. The situation may be different for relatively dense 
BEG, when the resulting ID model features deviations 
from the cubic interactions (for the binary condensate 
without the SO coupling this situation was elaborated in 
Ref. [HI), which may be a subject for separate analysis. 
The ID Hamiltonian gives rise to coupled GPEs, 
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+ (5i|^jf +5l03-,f J = 1,2, 



(1) 



potential /i correspond to ^{x,t) = e *'^*$(a;), where 
spinor 4> = col{0i, 02} obeys the system 
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(2) 



with 0j- = d(j)j/dx. In the free space (no trapping po- 
tential) and in the absence of the linear coupling, the 
commonly known condition for the immiscibility of the 
binary condensate is g > gi > [2^ [26|. The trap- 
ping potential and Zeeman splitting shift the transition 
to the miscibility from g = gi to larger values of g 271. 
Note that, unlike the equations considered in Ref. jl2| . 



model ID) is specific to the SO-coupled system, as it did 
not occur previously in fiber optics. In principle, there 
is a chance to implement this system in nonlinear op- 
tics too, but in a different context, using settings such as 
the recently reported Dirac model [l^ , which is obtained 
from Eq. Q by neglecting the kinetic-energy terms. 



III. ZERO ZEEMAN SPLITTING 

At = 0, Eqs. (H)) give rise to an evident 
solution, V'l = "02 = exp (iKx -t- iK^f/2) X (x, i), 
with function x obeying the standard GPE: idtX = 
[-(1/2)92 + + {gi + g) IxP] X- Below, we refer to 
this solution, characterized by nonzero superfluid veloc- 
ity, as a current state. At the same time, the sym- 
metry of system ([2]) with £7 = admits solutions of 
two other types: an odd mode, with 02(2^) = i0i(— a;), 
where 4>i{x) and 02 (a;) are purely real and imaginary 
functions, respectively; and an even mode, with 0i(x) 
and 02 (a:) having opposite parities, e.g., 0i(a;) is even 
and real, while 02 (a;) is odd and imaginary, see exam- 
ples in Fig. [H Odd and even modes correspond to sym- 
metric and anti-symmetric distributions of the pseudo- 
magnetization density, A4 = (|0i(a;)|2 — |02(a;)p) /N, 

whose integral M = A4dx ^ {Ni — N2)/N defines 
the total magnetization. 

Possible types of nonlinear modes, which constitute 
one-parameter families of solutions, can be identified by 
the analysis of bifurcations of the families from the lin- 
ear limit, — >■ (or, alternatively, gi = g = 0), which 
leads to the eigenvalue problem, Ti^ = jl^ (hereafter, 
the tilde stands for the linear limit). It is easy to see 
that the spectrum of this separable system is double- 
degenerate: fin = 71 + (1 — k2)/2, 71 = 0, 1, 2, ... Eigenvec- 
tors, = col{0i_„, 02, n}, can be chosen as arbitrary su- 
perpositions of two mutually orthogonal ((^•„.+ , ^n,-) — 



which conserve the total number of atoms, A^ = A^i + A^2, [2J]) spinors: ^„ = C'+*n.+ + C-^n,-, where $„,± 



with Ni^2 = |V'i,2r rfa;. Below we consider the 

generic situation {g ^ gi), and keep A^ as a free param- 
eter, fixing positive coefficients gi and g, which account 
for the intra- and inter-species repulsive interactions, re- 
spectively Stationary modes of Eq. ([T]) with chemical 



e=^"'^e~^"/2i7„(a;)col{l,±l}, and iJ„(a;) is the Hermite 
polynomial. 

However, an arbitrary set of constants C± in does 
not correspond to nonlinear eigenmodes. To determine 
the respective constraints, which implies lifting the degen- 
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FIG. 1: (Color online) Generic examples of nonlinear modes 
bifurcating from the lowest linear eigenstate with n — 0. Left 
[right] panels display odd modes for N = 8 [even modes with 
even <j}i{x) and odd (j>2{x) for A'" = 18]. Top [bottom] rows 
correspond to «: = 0.1 [k = 3]. Modes (b)-(d) are stable, while 
mode (a) is unstable. Other parameters are gi = g/2 — 1, 
Q = 0. Different vertical scales for (^i(x) and 4>2(x) in panel 
(b) indicate that |0i(2;)| ^ |<^2(a::)|. 



eracy by the repulsive nonlincarity, we introduce expan- 
sions ^j.^ fln + e'^fXn^ + . . . and * = e*„ + e^^n^ + . . . , 
which satisfy Eqs. © with ft = at the leading or- 
der with respect to e, which characterizes the strength of 
the nonlinearity. At order e^, one has — /i„)€>Jf'' = 



(2) A 



with 



F — 



(5l|01,nP 
(5l|02.nP 



' g\(t>l,n\'^)<t>2,n 



(3) 



(3) 

The solvability condition for results in two equa- 

tions, /i„ (*„,±,*„) = (*„,±,F,i), where constant 
C± are the unknowns. Focusing on the lowest linear 
eigenstate (the one which corresponds to n ^ 0), it 
is straightforward to find solutions of the latter equa- 
tions. One solution has C_ ~ and C+ 7^ 0, 
the respective linear mode being = 4'2,nix). 

In the nonlinear regime, it keeps the same structure, 
4>i{x) = 4>2{x), representing the current state, as de- 
fined above. Another solution has C+ = ±zC_. In 
the hnear limit, the associated eigenmode is $0 = 
g-K /2 {cos(7r/4 -I- Kx), i cos(7r/4 — ra)}, which extend 
into the above-mentioned odd nonlinear mode, with 
4>2{x) = iipii^x). Finally, the third solution, with 
C+ = ±C_ and ^0 = col {cos(Ka;), i sin(K2:)}, rep- 

resents even modes in the nonlinear regime. Importantly, 
no other solutions exists in the weakly nonlinear system. 

Thus we have identified three families of nonlinear 
modes bifurcating from the lowest eigenstate n = 0. Nu- 
merically found examples of the nonlinear modes from 
these families arc shown in Fig. [TJ While, as said above, 
one should expect the immiscibility due to the repulsive 
interactions, most modes feature the striped antiferro- 



magnetic structure, represented by alternating domains 
of states It) and \D , rather than two large domains sep- 
arated by a single (DW, which is the case in the conven- 
tional binary BEG [13, Hlj . It is seen that the number 
of the domains increases with the SO coupling strength, 
the single DW being present only in panel [Ija) . This 
patterning is precisely explained by the spatially peri- 
odic factors in the above analytical solutions. In quali- 
tative terms, the interplay of the HO trap with the SO 
coupling gives rise to scale / ~ 1/k that determines the 
periodicity of the striped patterns in Fig. [1] (the analyt- 
ical solutions yield / = tt/k). It is relevant to mention 
that multi-domain patterns were also reported in Ref. 
[2^ in the spinor (three-component) BEG under the ac- 
tion of external magnetic field (without the SO coupling). 
An essential difference of our system is that the transi- 
tion between different numbers of domains is controlled 
by the intrinsic strength of the SO coupling, rather than 
by an additionally introduced magnetic field. For the 
SO-coupled system with spin 2, two-dimensional multi- 
domain patterns were reported in Ref. [30| . but the sta- 
bility of those patterns (which is the main subject of the 
present analysis) was not studied. 

A similar analysis has been performed for higher-order 
nonlinear modes, i.e. for ones stemming from excited 
states of the linear system (n > 1, see above). They 
also feature the multi-DW structure of the odd and even 
types, with the number of domains increasing with the 
strength of the SO coupling, see the examples displayed 
in Fig. H 

In view of the coexistence of many multi-DW patterns, 
their stability is a crucial issue. Here we are meaning the 
experimentally relevant dynamical stability, determined 
by the Bogoliubov-de Gennes (BdG) equations derived 
from Eqs. ([1]), rather than the thermodynamic stabil- 
ity. The spectrum of the BdG equations was found by 
means of standard methods for the solution of the cor- 
responding eigenvalue problem [26j , and the results were 
verified by direct simulations of the perturbed evolution 
of the modes. A surprising conclusion is that the sim- 
plest pattern with the single DW, shown in Fig. [Ija), 
is unstable^ on the contrary to the situation in the ordi- 
nary binary BEG [l^,!!!], while the multi-DW patterns, 
displayed in panels (b-d), are stable. This finding may 
be explained by the fact that stable are those structures 
which comply with the above-mentioned spatial scale I 
naturally selected by the system. In direct simulations of 
the perturbed evolution, see Fig. (3]), the unstable single- 
DW state spontaneously develops strong oscillations, in 
agreement with the oscillatory instability predicted by 
the BdG analysis. The resulting formation of robust dy- 
namical modes in the form of breathers is a physically 
relevant result too. 

Furthermore, we have found that higher-order nonlin- 
ear modes, with the number of domains increasing with 
strength of the SO coupling, which bifurcate from the 
excited linear states, with /i = /ii,2,... (see above), may 
also be stable in the nonlinear system, see examples of 
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FIG. 2: (Color online) Examples of stable nonlinear modes 
bifurcating from the first excited {n = 1) linear eigenstate 
(rather than from the ground state with n — 0, cf. Fig. [T] 
Panels (a) and (b) show, severally, odd and even modes for 
51 =3/2 = 1 and ^ = 0. 




FIG. 3: (Color online) The evolution of the unstable nonlinear 
mode shown in Fig.[T](a). Under the action of the instability, 
the mode spontaneously transforms into a persistent oscilla- 
tory state. 
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FIG. 4: (Color online) Stability domains (shaded) for odd 
(a) and even (b) nonlinear modes bifurcating from the lowest 
linear eigenstate, with n = 0, in the system with the SO 
coupling in the absence of the Zeeman splitting, (c) and (d) ; 
The overlap integral, K, versus the total norm. A'', for odd 
(c) and even (d) modes. Curves labeled -ffo.i.i.s correspond 
to K = 0.1, 1, and 3, respectively. The dashed curve in (d) 
shows the total magnetic moment, M, as a function of A'^ for 
K = 1. For K = 0.1 and k = 3, M does not vary significantly 
with A'^. Other parameters are gi = g/2 = 1, Q — 0. 



stable higher-order modes in Fig. [2j 

The full stability chart for the odd and even nonlin- 
ear modes originating from fj, = fiQ is fairly complex, as 
shown in Fig. HI For sufficiently small values of k and 
in the linear limit iV — >■ only the even modes are sta- 
ble. This is readily explained by the fact that, at small 
K, the system admits an obvious even configuration with 
almost all atoms falling into a single state, see Fig. [Ijb) . 
For larger k, the stability diagram features a zebra-like 
structure, with alternating stability regions of the even 
and odd modes (while the aforementioned current states 
are completely unstable, cf. Ref. [3l|). The stability ar- 
eas for the odd and even modes approximately (but not 
precisely) complement each other, which is a consequence 
of the competition between coexisting nonlinear modes 

These stability results arc the most essential findings 
reported in this work, as they demonstrate the previ- 
ously unreported transition to complex ground states in 
the system, under the action of the SO coupling. In fact, 
this transition may be expected in diverse systems be- 
yond the model of the linearly-coupled binary BEG. 

The stable patterns are naturally characterized by the 
total pseudo- magnetization M {M = ±1 means that the 
condensate is in a single-domain state), and by the mis- 



cibility factor, 

/ + 00 / P + CO 

l^pil^lM^dx / / (|0ip + |02p)'dx, (4) 

which takes values < A' < 1, with K = and K = 1 
corresponding, severally, to the completely immiscible or 
miscible state. Note that the immiscibility, K 0, may 
be achieved not only if the phases are spatially separated, 
but also if one phase disappears, and the condensate falls 
into a single-domain state, as in the case of the even 
mode in Fig. [1] (b). In the linear limit N 0, one has 
M = e~'^ for the even eigenstates, and M = for the 
current and odd ones, while AT = (1 =F e^^'^ )/2 for the 
odd (— ) and even (-I-) states. 

Figure ^c) shows that the nonlincarity enhances the 
immiscibility of the odd modes, although the effect may 
be significantly reduced by the coupling, as curve 
demonstrates. The effect of the interplay between the 
nonlincarity and SO coupling is more sophisticated for 
even modes, see Fig. Hl^d). When the SO coupling is 
sufficiently weak (curve ifo.i) or strong (curve K^), the 
nonlinear modes preserve properties of their linear coun- 
terparts: at K = 0.1 the mode is nearly fully polarized 
(A/q.i ~ 1, Kq.i ^ 1), while the mode with k = 3 
(A/3 ^ 1, K3 ?a 0.5) represents a partially miscible state. 
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where both components are equally represented, with Ai 
close to zero. 

The most interesting behavior is observed at interme- 
diate values of the SO coupling, which corresponds to 
K = 1 in panel (d). In this case, by changing norm 
N one draws the condensate into the immiscible state 
(curve i^i), and, at the same time, changes M between 
zero and its extreme value +1. Bearing in mind that the 
BEG nonlinearity may be controlled via the Feshbach res- 
onance [2^, the observed behavior suggests a possibility 
of switching between the different phases with the help 
of external fields. 



IV. SYMMETRY BREAKING DUE TO THE 
ZEEMAN SPLITTING 

In the linear limit N ^ the Zeeman splitting [f2 7^ 
in Eqs. ([T])] removes the energy degeneracy. Now, the 
system does not admit current states and odd modes, 
while stable even modes, characterized by even 0i(a;) and 
odd 02 (a^) (or vice versa), persist at f2 7^ 0. The evolu- 
tion of the eigenstates under the action of the increasing 
Zeeman splitting, il, is illustrated by the bifurcation di- 
agram in the plane of (Af, fi), shown in Fig. [5j Starting 
with the point where an odd (even) mode is stable (unsta- 
ble) at J7 = 0, we observe a transformation of the stable 
odd mode into stable asymmetric ones [see an example 
in Fig. [S] (a)], whose branch approaches the branches of 
the even modes. The instability of the even modes per- 
sists up to the point where the branches of even and 
asymmetric modes merge (points Pi,2)- At these points, 
the asymmetric modes disappear and the symmetric even 
ones become stable [an example is shown in Fig. EKb)]. 
Thus, a pitchfork bifurcation occurs at points -Pi, 2, which 
is a typical example of the spontaneous symmetry break- 
ing (or restoration) [s^ . 

The bifurcation diagram in Fig. [5] suggests that, by 
varying the Zeeman field which induces the Zeeman 
splitting, one can perform a controllable switch between 
two opposite magnetizations in the SO-coupled binary 
system. This suggestion is confirmed by simulations, as 
shown in Fig. [51 




6 



FIG. 5: (Color online) The bifurcation diagram (the left 
panel) for «; = 1 and A'' = 5, in the presence of the Zeeman 
splitting, Jl, the other parameters being as in Fig. [l] Stable 
(unstable) modes correspond to solid (broken) fragments of 
the curves. Panels (a) and (b) show profiles of stable modes 
marked by points Pa,b. 




X X 

FIG. 6: (Color online) The evolution of a nonlinear mode 
subjected to an adiabatic change of Q.. Panel (a) shows the 
input and output spatial profiles. Panels (b) and (c) display 
intensity plots. The input at f = is a stable odd mode 
at n = 0, the other parameters being as in Fig. [5] In the 
course of the simulation, Q. was adiabatically changed from 
f7 = Oatt = Otor2 = -0.25 at t = 10^. As a result, the 
condensate switches from the initial odd mode to an even one, 
as predicted by the diagram in Fig. O 



V. CONCLUSION 

We have reported the existence of even, odd, and asym- 
metric nonlinear modes in the effectively ID self-repulsive 
binary BEG with the SO and Zeeman splitting, confined 
by the axial HO potential. The interplay between the 
coupling and the potential gives rise to the modes featur- 
ing alternating domains with opposite directions of the 
pseudo-spin, which is explained analytically in the case 
of the weak nonlinearity. Noteworthy findings are the 
stability of the multi-DW patterns, while the one with 
the single DW is unstable, and the stability alternation 
between the even and odd structures. The implication 



of these results is the transition from simple to complex 
ground states, driven by the SO coupling. The inclu- 
sion of the Zeeman splitting results in the transformation 
of even modes into asymmetric ones, which suggests a 
possibility of controllable switching between states with 
opposite pseudo-magnetization. These effects, including 
the transition to the complex ground states, may be as 
well expected in other physical setting emulated by the 
SO-coupled binary condensates, such as solid-state media 
and bimodal guided-wave propagation in optics. 
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